Structure, rotational dynamics, and superfluidity of small OCS-doped He clusters 
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The structural and dynamical properties of OCS molecules solvated in Helium clusters are studied using 
reptation quantum Monte Carlo, for cluster sizes n = 3 — 20 He atoms. Computer simulations allow us to 
establish a relation between the rotational spectrum of the solvated molecule and the structure of the He solvent, 
and of both with the onset of superfluidity. Our results agree with a recent spectroscopic study of this system, 
and provide a more complex and detailed microscopic picture of this system than inferred from experiments. 
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Solvation of atoms and molecules in He nanodroplets pro- 
vides a way to study their properties in an ultra-cold matrix, 
and also offers a unique opportunity to probe the physics of 
quantum fluids in confined geometries. Research in this field 
has been recently reviewed in Ref. with emphasis on ex- 
periment and hydrodynamic modeling, and in Ref. [^], with 
emphasis on computer simulations. 

Carbonyl sulfide (OCS) is one of the most widely studied 
dopants of He clusters (OCS@He„) ^, both because of 
its strong optical activity in the infrared and microwave spec- 
tral regions, and also because the OCS@Hei complex is spec- 
troscopically well characterized [H ||], thus providing a solid 
benchmark for the atom-molecule interaction potential which 
is the key ingredient of any further theoretical investigation. 

In the quest of fingerprints of superfluidity in the spectra 
of small- and intermediate-size He clusters. Tang et aJ. have 
recently determined the vibrational and rotational spectra of 
OCS@He„ at high resolution for n = 2 - 8 [||]. The main 
results of that investigation are: j) the rotational constant of 
the solvated molecule roughly equals the nanodroplet (large- 
n) Umit at n = 5, but then undershoots this asymptotic 
value up to the maximum cluster size (n — 8) attained in that 
work; ii) the centrifugal distortion constant has a minimum at 
n = 5, thus indicating that the complex is more rigid at this 
size; in) the fundamental vibrational frequency of OCS is not 
a monotonic function of the number of He atoms, but it dis- 
plays a maximum at n = 5, again suggesting a stronger rigid- 
ity at this size. Findings ii) and Hi) suggest — and our study 
of the rotational spectra confirms — that n = 5 is a magic size 
related to the structure of solvent atoms around the solvated 
molecule. Finding i) implies the existence of a (yet to be de- 
termined) minimum in the rotational constant as a function of 
the cluster size. The occurrence of this minimum was inter- 
preted as due to quantum exchanges which would decrease the 
effective inertia of the first solvation shell and would thus be a 
signature of the onset of superfluidity in this finite system [|^]. 

Quantum simulations are complementary to the experiment 
for understanding the properties of matter at the atomic scale. 
In fact, while being limited by our incomplete knowledge of 



the inter-atomic interactions and by the size of the systems 
one can afford to examine, simulations provide a wealth of de- 
tailed information which cannot be obtained in the laboratory. 
In this paper we show how state-of-the art quantum Monte 
Carlo simulations can be used to get an unparalleled insight 
into the exotic properties of small He clusters. At the same 
time, some of the dynamical properties thus predicted are so 
sensitive to the details of the atom-molecule interactions, to 
provide a very accurate test of the available model potentials. 

The structure and the rotational dynamics of OCS@He„ 
clusters are studied here in the small-to-intermediate size 
regime (n = 3 — 20), using reptation quantum Monte Carlo 
(RQMC) a technique by which ground-state properties 
(such as the density distribution and various static and dy- 
namic correlations) can be determined with high accuracy. 
In particular, the analysis of the dynamical dipole correlation 
function allows us to resolve, for the first time theoretically, 
several rotational components of the spectrum of an interact- 
ing quantum system. RQMC combines the best features of 
the path-integral (PIMC) and diffusion quantum Monte Carlo 
(DQMC) techniques, yet avoiding some of their drawbacks: 
RQMC is not affected by the systematic errors which plague 
the estimate of ground-state expectation values in DQMC (the 
so called mixed-estimate and population-control biases) and 
gives easy access to time coiTelations, from which dynam- 
ical properties can be extracted [^; at the same time, the 
RQMC technique is specially tailored to the zero-temperature 
(ground-state) regime, where PIMC becomes inefficient. 

Both the He-He and the He-OCS interactions used here are 
derived from accurate quantum-chemical calculations |^, [lo| ]. 
The OCS molecule is allowed to perform translational and 
rotational motions, but it is assumed to be rigid. The trial 
wavefunction is chosen to be of the Jastrow form: <I>t = 



exp 



where 



is the dis- 



i<j 

tance between the i-th and the j-th helium atoms, ri = |rj| 
the distance of the i-th atom from the center of mass of the 
molecule, and 6i is the angle between the molecular axis and 
Ti. On account of the high anisotropy of the He-OCS poten- 
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FIG. 1: Logarithm of the multipole correlation function, norm- 
lized by the spectral weight of the lowest-lying state: e(r) — 
— log {cj{t)/Aj) /t, for 71 = 5 (lower panel) and n = 10 (upper 
panel). Time units are °K^^. 

tial, Ui is expanded as a sum of five products of radial func- 
tions times Legendre polynomials. All radial functions (in- 
cluding U2) are optimized independently for each cluster size 
with respect to a total of 27 variational parameters. The prop- 
agation time is set to 1 °K^^, with a time step of 10"'^ °K^^. 
The effects of the length of the time step and of the projection 
time have been estimated by test simulations performed by 
halving the former or doubling the latter. These effects were 
barely detectable on the total energy, and negligible with re- 
spect to the statistical noise for the properties on which we 
base our discussion below. 

The absorption spectrum of a molecule solvated in a 
non polar environment is given by the Fourier transform 
of the autocorrelation function of its dipole, d: I{uj) cx 
2^E„ \{^o\d\^n)\^6iEn-Eo~Lu) = /e^-*(d(i).d(0))di, 
where '^q and "ifn are ground- and excited-state wavefunc- 
tions of the system respectively, Eq and En the correspond- 
ing energies, and (.) indicates ground-state expectation val- 
ues. The dipole of a linear molecule — such as OCS — is ori- 
ented along its axis, so that the optical activity is essentially 
determined by the autocorrelation function of the molecular 
orientation versor: c{t) = {n{t) ■ n(0)). Within RQMC 
the analytical continuation to imaginary time of this corre- 
lation function, c(t) — c{-~iT), can be conveniently calcu- 
lated without any other approximations than those implicit 
in the use of a given parametrization of the interatomic po- 
tentials From now on, when referring to time correla- 
tion functions, we will mean quantum correlations in imag- 
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FIG. 2: Rotational constants (upper panel) and rotational distortion 
constant (lower panel, see text) of OCS@He,i, as functions of the 
cluster size, n. Triangles indicate the results of the present simula- 
tion, whereas circles are experimental data from Ref. |b|]. 



inary time, and will we drop the bar over c(t). Analyt- 
ical continuation to imaginary time transforms the oscilla- 
tory behavior of the real-time correlation function — which is 
responsible for the (5-like peaks in its Fourier transform — 
into a sum of decaying exponentials whose decay constants 
are the excitation energies. En ~ Eq, and whose spectral 
weights are proportional to the absorption oscillator strengths, 
(4'o|d|^'„) p. Dipole selection rules imply that only states 
with J — 1 can be optically excited from the ground state. 
Information on excited states with different angular mo- 
menta, J, can be easily extracted from the multipole corre- 
lation functions, cj(t), defined as the time correlations of 
the Legendre polynomials: cj(t) = (P7(n(T) • n(0))) = 

( 2JTT Em=-j YJm {n{r))YjM (n(0))). For each value of 
J, the value of the lowest-lying excitation energy, e,j — i.e. 
the smallest decay constant in cj(r) — as well as the corre- 
sponding spectral weight, Aj, can be extracted from a fit to 
c,/(t). We have verified that the values so obtained are rather 
insensitive to the details of the fitting procedure. 

In Fig. 1 we display multipole correlations calculated up 
to J = 5 for OCS@He5 and OCS@Heio. For a rigid-top, 
time correlations are described by a single exponential for 
each value of J, cj(r) = e"'^-''^, whose decay constant is 
ej = BJ{J + 1). The rigid-top energy levels are displayed 
as dashed lines. Time correlations deviate from a single- 
exponential behavior only at short times, thus indicating that 
the non rigidity of the cluster only affects the spectra at very 
high frequencies, the deviation being larger the larger the an- 
gular momentum. The dependence of e on J is accurately 
predicted by the rigid-top model for n = 5, while it is less 
so for n = 10. In all cases the excitation energies are well 
described by a formula containing an effective centrifugal dis- 
tortion constant: ej = BJ{J + 1) ~ D.P{J + 1)^ Q. 

In Fig. 2 we display the rotational constant, B, and cen- 
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trifugal distortion constant, D, as functions of the cluster size. 
B displays a minimum at n = 8 — 9, while D is mini- 
mum at n = 5, indicating a greatest rigidity of the He@OCS 
complex at this size. In order to discuss these findings and 
their relation with the structure of the cluster, with the on- 
set of superfluidity, and with the quality of the interaction 
potential used for the simulation, we first define p(r) as the 
He number-density distribution, and 0(r) as the ground-state 
expectation value of molecule-atom angular motion corre- 
lation: (p{r) ~ — (^o|L ■ l(r)|5'o). In this expression L 
indicates the angular momentum of the OCS molecule and 
l(r) is the angular-current operator of He atoms at point r: 

l(r) = -^nr xj:i{i:Sir-ri) + S{r-ri)£^). The 

angular correlation function, (/)(r), is non-negative and it is 
large whenever the angular motion of He atoms at point r is 
strongly correlated with the rotation of the solvated molecule, 
thus contributing to the molecular effective moment of inertia. 

In Fig. 3 we display p(r) and (j>{r), as calculated for clus- 
ter sizes n — 5,8, 10, 15. The He density reaches its maxi- 
mum where the He-OCS interaction potential is minimum, i.e. 
in a doughnut surrounding the molecule perpendicular to the 
OC bond. This doughnut can accommodate up to 5 He atoms 
which rotate rather rigidly with the OCS molecule, as demon- 
strated by the large value of the angular correlation function. 
Each He atom within the doughnut gives a same contribution 
to the cluster moment of inertia, resulting in a constant slope 
of B vs. n for n < 5. 

Secondary minima of the He-OCS potential exist on the 
molecular axis near the two poles, the one nearest to the S 
atom being deeper. For n > 5 He atoms spill out of the re- 
gion of the main minimum and spread towards other regions 
of space, preferentially near the molecular poles. A detailed 
analysis of the He density plots for n — 6, 7, 8 shows that — 
due to quantum tunneling to and from the doughnut — ^both 
polar regions start to be populated as soon as the doughnut oc- 
cupation is completed (i.e. forn > 5). The number of Helium 
atoms near the sulfur pole is 0.2, 0.3, and 0.7, for n = 6, 7, 8, 
and 0.5, 0.8, and 1.1 near the oxygen pole. Helium density 
is larger near the oxygen pole — ^rather that near sulfur where 
the potential is more attractive — because a smaller energy bar- 
rier and a smaller distance from the absolute minimum in the 
doughnut make quantum tunneling between the absolute and 
the secondary minima easier in this case. 

Quantum tunneling is the key for understanding the onset of 
superfluidity, as well as the sensitivity of rotational spectra to 
the details of the He-OCS interaction. For n = 6, 7, 8 a size- 
able He density is found not only near the molecular poles, 
but also in the angular region between the doughnut and the 
oxygen pole, where the energy barrier is small. Non negligi- 
ble angular correlations indicate that atoms in this region (and 
not only near the poles) contribute to the effective moment of 
inertia. For n > 8 He atoms start to fill the angular region 
between the doughnut and the S pole (see the n = 10 panel of 
Fig. 3). The closure of He rings in the sagittal planes makes 
atomic exchanges along these rings possible, thus triggering 



the same mechanism which gives rise to superfluidity in in- 
finite systems [|ll]]. The onset of superfluidity is responsible 
not only for negligible angular correlations in this region, but 
also for the decrease of these correlations in the doughnut re- 
gion where the largest contributions to the cluster inertia come 
from. This is demonstrated in the bottom-right panel of Fig. 
3, where one sees that the maximum of the angular correlation 
in the doughnut region starts decreasing for n > 8. 

It is also interesting to examine the dependence upon clus- 
ter size of the spectral weights, Aj. It turns out that, while for 
J = 1 the spectral weight displays a maximum at n = 4 — 5 
and a minimum at n = 8, the maximum at n = 4 — 5 is not 
followed by any minimum for J = 5. We interpret the max- 
imum of Aj at n — 4 — 5 as another manifestation of the 
greatest rigidity of the He@OCS cluster at this size, while the 
behavior for n > 5 is determined by the number of solvent 
states available to couple with the lowest-lying molecular ro- 
tation. This number increases with the cluster size for J — 5, 
while it has a maximum at n = 8 — 9 for J = 1 . These find- 
ings are consistent with the emergence of superfluid behavior 
at cluster size, as characterized by the reduction of low-lying 
excited states in the He matrix rnTIl . 



This quantum tunneling makes the rotational constants for 
n > 5 rather sensitive to the height and width of the barriers 
between different minima of the OCS-He interaction poten- 
tial. On the other hand, the details of the potential off the 
main minimum in the doughnut hardly affect the spectra for 
n < 5 and, in particular, that of OCS@Hei which is used 
as a benchmark for the quality of the potential [10|. Our 
simulations predict rotational constants which are in excel- 
lent agreement with experiment for n < 5, while for n > 6 
the predicted values of B are too large and display too small 
a slope as a function of n. We believe that both these facts 
are due to the smallness of the barriers between the main and 
the secondary energy minima, which determines an excessive 
spill-out of He atoms off the main minimum for n = 6 and 
an excessive propensity towards atomic exchanges for n > 6. 
We expect that larger energy barriers would result in a better 
agreement between theory and experiment for 6 < n < 8 
and, possibly, in a shift of the predicted minimum of the rota- 
tional constant towards larger sizes. As a qualitative test, we 
have redone some of our simulations with a uniformly scaled 
potential (Vhc-ocs ^ aVno-ocs)- Although not very real- 
istic, this transformation captures the essential modifications 
to be done to the potential, i.e. it modifies the height of the 
barriers while leaving the location of the minima unchanged. 
For a = 1.2, we find that the values of B for n < 5 are left 
practically unchanged, while the agreement with experiment 
is much improved for n = 6, 7, 8. The resulting minimum 
value of B stays located near n = 8 — 9. This lets us believe 
that our prediction of the position of the minimum rotational 
constant as a function of the cluster size should be rather ro- 
bust with respect to possible improvements of the He-OCS 
effective potential. 



4 



FIG. 3: Upper left panel: He-OCS interaction potential. The atoms constituting the molecule are displayed as circles whose radius is the 
corresponding Van der Waals radius (sulphur on the left). Panels with red background: upper mid panel, He number density; lower mid panel, 
angular motion correlation function (see text). Colors go from red to violet as in a rainbow, as the value of the function increases. Lower right 
panel: maximum of the angular motion correlation function, as a function of the cluster size. 




In conclusion, we have shown that state-of-the art quantum 
Monte Carlo simulation provide both a qualitative explanation 
of the relation existing among structure, dynamics, and super- 
fluidity in small He clusters, and a sensitive test of the quality 
of atom-molecule potentials, in regions which hardly affect 
the spectra of the monoatomic molecular complex. 
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